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ABSTRACT 

The operation of replenishment at sea is investigated for three 
supply ships and L combatants using queueing theory concepts and a 
randan walk model in three dimensions. The distribution for total 
replenishment time given the initial number of canbatants to be 
replenished by each supply ship and a specified cyclic order of 
replenishment is expressed in terms of its Laplace transform under 
the assumption of independent, exponential service times for each 
supply vessel. A method for counting all possible sequences of 
replenishment is not found, but sane preliminary counting techniques 
are developed which may be useful in its eventual determination. 
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I. INTRODUCTION 



A major aspect of modern naval policy is the stationing of 
task forces throughout the world, and the supply and replenishment 
of these task forces is a critical problem with which naval plan- 
ners are faced. It is vital for ships at sea to be provided with 
a reliable flow of supplies for then to maintain a continuous 
readiness posture. 

Forces afloat are presently supplied with fuel, food, ammunition, 
and other bulk goods primarily by large supply ships which accompany 
task forces during their operations and may periodically proceed to 
bases ashore to replenish themselves . Transfer of goods at sea is 
largely accomplished by combatants coming alongside these support 
ships and receiving supplies by use of various rigs. The nature of 
the primary alongside method of transfer and its consequent grouping 
of numbers of ships over a period of time increases the vulnerability 
of all participating vessels. Therefore, time of replenishment of 
the whole combatant task force is of primary concern. The problem 
of scheduling a replenishment operation of a combatant task force 
by a group of supply ships such as oilers, ammunition ships, and 
other support vessels has been studied by McCullough [3] , Gordon 
and Copes [2] , Patterson [5] , and Waggoner [6] . Waggoner built a 
mathematical model to study the scheduling of a fixed number of 
combatants by two supply ships and computed the distribution of total 
replenishment time of all combatants. It is the purpose of this paper 
to consider the analogous problem with three supply ships. 
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II. FORMULATION OF THE PROBLEM 



In general the problem may be considered frcm a queueing theory 
approach with each supply ship representing a service facility and 
each combatant a cus toner. It may or may not be necessary for each 
customer to visit each facility, depending upon the needs of the 
customer and the number of duplicate facilities. Waggoner [6] and 
Milch and Waggoner [4] employ a random walk model for study of a 
two facility queue which operates in parallel and series simultaneously 
with a finite number of customers. Such a model is applicable to the 
replenishment situation. It is the purpose of this paper to extend 
the random walk model to the case of three supply ships under the 
following assumptions: 

(i) A finite number of combatants are initially present, 

and each must be served once by all three supply ships. 

(ii) The service times at each supply ship are independent, 
exponential random variables with parameters A, y, and 
v for ships I, II, and III, respectively, and each one 
serves only one combatant at a time. 

(iii) Service is simultaneous and consecutive. When a combatant 
has been replenished in one queue he moves to the next 
queue in cyclic order. Each facility continues in operation 
until it has served all combatants. 

(iv) No additional combatants join the queues from outside 
the system after the commencement of operations. 

(v) When a combatant has been served by all three supply 
ships he departs the system. 
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(vi) Transit times between queues are considered negligible. 
Figure 1 illustrates the initial positions of all vessels and 
the order of flow fran one queue to another. 




M + N + P = L 
Figure 1 
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III. THE THREE DIMENSIONAL MODEL 



Milch and Waggoner [4] have represented an analagous two facility 
queueing systen by a randan walk in a two dimensional lattice. The 
random walk begins at the origin at time zero and proceeds by horizontal 
and vertical steps to the point (L,L) where L = M + N, and M and N are 
the initial number of cus toners in queues I and II, respectively. A 
horizontal step corresponds to completion of servicing in one queue, 
and a vertical step is a canpleted service in the other. Hawevera 
all paths to (L,L) are not possible since if X(t) and Y(t) represent 
the amount of custaners served in each respective queue at time t, 
then X(t) < M + Y(t) and Y(t) < N + X(t) . Figure 2 shews a path in 
such a representation. 




Figure 2 

In the two dimensional model the distribution of total operation 
time, T, is found in terms of its Laplace transform. 
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If three supply ships are present as in Figure 1 a three 
dimensional random walk may be formulated which begins at the 
origin at time zero and proceeds by unit steps at a time in one 
of the three axial directions to the point (L,L,L) where L = M + N + P, 
and M, N, P are the initial number of ships in the first, second, and 
third queue, respectively. All paths to (L,L,L) do not represent a 
replenishment operation, however, since if X(t) , Y(t), and Z (t) are 
the number of ships served in queues I, II, and III respectively, 
then X(t) < M + Z (t) , Y(t) < N + X(t) , and Z (t) <. P + Y(t) are 



constraints that must hold during the whole operation. These con- 
straints define boundary planes for the three dimensional lattice. 

The random walk paths are further bounded by X(t) = L, Y(t) <= L, 
Z (t) = L. For the purpose of clarity in further development these 
planes and their intersections will be catagorized as follows: 



Boundary 

A 

B 

C 

D 

E 

F 



Constituents 



Note: 



x = z+ M, z = L, and their intersection 
y=x+N, y=L, and their intersection 
z=y+P, z=L, and their intersection 
Intersection of A and B 
Intersection of A and C 
Intersection of B and C 
D is not regarded as part of A or B. Similar statements may 
be made about E and F. 
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IV. DISTRIBUTION OF TOTAL REPLENISHMENT TIME 



A. PATH PROBABILITIES 

The first move of the random walk is made one unit in an axial 
direction with probabilities A/A + y + \>, y/A + y + v, and v/A + y + v 
for axes X, Y, and Z, respectively, since these are the probabilities 
of stations I, II, or III, respectively, completing replenishment of 
the first combatant under the assumption of independent, exponential 
service times. Due to the memoriless property of the exponential 
distribution these probabilities are true at all steps along the random 
walk, except at points on the boundaries whose probabilities are as 
follows : 



Boundary 


Axial Direction 


Probability 


A 


Y 


y/y + v 




Z 


v/y + v 


B 


X 


A/A + v 




Z 


v/A + v 


C 


X 


A/A + y 




Y 


y/A + y 


D 


Z 


one 


E 


Y 


one 


F 


X 


one 



The total time for the replenishment is the time needed for the 
random walk to proceed from the origin to the point (L,L,L) which is 
the sum of the times needed to make each of the 3 L steps required for 
the walk. The distribution of time for any of the steps not originating 
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on a boundary is the minimum of three independent, exponential 
distributions and therefore is itself exponential with parameter 
X + y + v . Similarly the duration of steps originating on bound- 
aries A, B, and C is exponential with parameters y + v , X + v, and 
X + y , respectively . The distribution of steps originating on bound- 
aries D, E, and F is exponential with parameters v, y, and X 
respectively . 

Let Q(i, j ,k,£,m,n) be a route that has exactly i, j, k. £, m, n 
positions on boundaries A, B, C, D, E, F, respectively. Then 
Q(i, j ,k,£,m,n) has 3L - (i + j+ k + £ + m + n) positions along the 
path that are not on a boundary. The total replenishment time, T, 
is composed of the following random variables, given that Q(i, j ,k,£,m,n) 
occurs: 



T = A. + B.+C. + D. + E + F + G (4-1) 

i j x £ m n g 

where A^ is the sum of i steps each of which originates on boundary A, 

and B., C, , D„ , E , and F are similarly defined. G is the sum of 
j k' £ m n J g 

g steps none of which originates on any boundary where g = 3 L - i - 
j-k-fc-m-n. 

Path Q(i, j ,k,£,m,n) contains L steps in each of the three axial 
directions. If the L steps in the x direction are examined, exactly 
n of them will begin on boundary F each with probability one. Seme 
number, j^, will begin on boundary B each with probability X/X + v 
and seme number, k^, begin on boundary C each with probability X/X + y. 
The remaining L-n-j^-k^ steps will begin from a non-boundary 
position each with probability X/X + y + v . Similar statements may be 
made concerning the steps in the other two axial directions, and the 
results are tabulated below. 
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Axis 


Number 


Step 


Direction 


of Steps 


Probability 


X 


n 


one 




A 


A/A + v 




1 — i 


X/X + y 




L - (n + j 1 + kj.) 


X/X + y + \) 


y 


m 


one 




i i 


y/y + v 




k 2 


y/A + y 




L - (m + i^ + k 2 ) 


y/A + y + v 


z 


£ 


one 




1 2 


v/y + v 




A 


\>/X + v 




L - (£ + i 2 + j 2 ) 


v/A + y + v 


where i = 


A + A' j = j l + j 2 ' k 


= Then 



P (Q) = r(i / j,k,i,m,n) ( y. A;) L (nf A +k l ; (tt -) 31 



, A , k l , p , L - (m+j-,+k 9 ) 

l A+ M ; ^ A+p+v 



(-H- -) jl 
\+\> J 




k 

2 



( v \ L - ( ^+i ? +j o ) / v \ 1 2 

A+p+v 2 p+v 



(— ) J 2 

v A+v 



P(Q) 



r(i / j,k,£,m,n) 



. A+p+v , i , A+p+v , j 
' p+v ' ' A+V 1 



. A+p+v v k . v ■ L-£ , p , L-m , A , L-n 
A+p A+p+v A+p+v A+p+v 



( 4 - 2 ) 



where Q = Q(i, j ,k,£,m,n) and r (i, j ^,£,111,^ is the number of paths 
from the origin to the point (L,L,L) which have i,j,k / £,m,n points 
in ccrrmon with boundaries A, B, C, D, E, F, respectively. 
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B. TRANSFORM OF THE DISTRIBUTION 



Frcm equation (4-1) the following Laplace transform can be 
tabulated : 

Randan Variable Laplace Transform 



A. 

1 



(p+v/p+v+s) 1 



B. 

3 






D 



£ 



E 



m 



G 

g 



(A+v/A+v+s) ^ 
(A+p/A+p+s)^ 
(v/v+s) ^ 

(p/p+s) m 

(A/A+s) n 

(A+p+v/A+p+v+s)^ 



The Laplace transform of the total replenishment time, T, is: 



f* T (s) = E(e ST ) = I E(e~ ST |Q) p(Q) 

i,j,k 

£,m,n 



* 

f (s) = E r(i,j,k,£,m,n) 
i/j/k, 

£,m,n 



, Apv 1 ,A+p+v 3L 
A+p+v' A+p+v+s^ 



, A+p+v+S' i , A+p+v+s j . A+p+v+Sx k 
p+v+S A+v+s ' A+p+s 



,A+p+v+s,. £ ,A+p+v+Ssm , A+p+v+s^ n 
‘ v+s ' ' p+s J ( A+s ' 
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v - COUNTING TECHNIQUES 



A. DISCUSSION 

An expression for r (i, j ,k,£,m,n) has not been developed in this 
paper. Hcwever, this chapter demonstrates several counting techniques 
for the two supply ship, two dimensional random walk which may be use- 
ful in the eventual determination of r(i, j ,k,£,m,n) . Additionally, 
for the three dimensional model an expression is found for the number 
of paths proceeding fran the origin to (u,v,w) which touch no bound- 
aries. In what follows, the binomial coefficient, (£) = n!/rl (n-r) i 

B. TWD DIMENSIONAL COUNTING 

1. The Number of Paths Not Touching Boundaries 

First a method will be formulated for counting the number of 
paths proceeding from (0,0) to (u,v) which do not touch or cross either 
boundary y = x + a, (a > 0) , or boundary y = x + b, (b < 0) . This 
number is equal to the total ( U + V ) paths going to (u,v) less the number 
of paths which touch or cross either boundary. By the Reflection 
Principle (see page 70 of [1] ) , the number touching or crossing y=x + a 
is and the number touching or crossing y = x + b is (y-b) * S: '- nce 

in both of these quantities paths touching or crossing both boundaries 
are included, the number of such paths must be added back. This quantity 
may be ccmputed by dividing it into two mutually exclusive categories — 
those paths which touch or cross boundary y = x + a first, and those 
paths which touch or cross boundary y = x + b first. Such paths will 
be counted by a double application of the Reflection Principle. In 
Figure 3 consider path 1 which falls into the first of the two categories. 
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The first reflection is with respect to line y = x + a. Both the 
boundary y = x + b and that part of path 1 beyond its first ccmrton 
point with y = x + a are reflected. The reflected path, path 1 ' , 
ends in point (v-a, u+a) . The common points of path 1 and y = x + b 
are reflected into ccrrmon points between path 1 ' and the line 
y = x + 2a-b (this being the reflection of line y = x + b) . Therefore, 
the number of paths in the first category is equal to the number of 
paths ending in (v-a, u+a) and touching or crossing the line y = x + 2a-b. 
The latter number is 



v-a+u+a . _ ,u+v . 
u+a- ( 2a-b) ^ u-a+b 



(5-1) 



by a second application of the Reflection Principle. Similarly the 
number of paths going to (u,v) after touching y = x + b first is 
^ u+a-b ^ * Suc ^ a i s path 2 in Figure 3 . Therefore the number 

of paths to (u,v) touching neither y = x + a nor y = x + b is 



,U+Vv _ ,U+V) 

{ u ; l u+a 



^ “ © + © + 



u+a-b 



u-a+b 



(5-2) 



By extending the boundary reflection employed in the preceding 
paragraph it is possible to compute the number of paths proceeding to 
(u,v) which touch or cross y = x + a, a > 0, first, y = x + b, b < 0, 
second, and y = x + c, c > a, third, before going to (u,v) . In Figure 
4 since y = x + a is the first boundary to be struck, reflect y = x + b 
and y = x + c with respect to y = x + a yielding y = x + 2a-b ard 
y = x + 2a-c, respectively. Since y = x + 2a-b is the next boundary 
to be struck by the reflected path, reflect y = x + 2a-c with respect 
to y = x + 2a-b, yielding y = x+2a - 2b+c. The twice reflected path 
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Figure 3 
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(v-a+b-c, u+a-b+c) 




Figure 4 



IS 



will then strike y - x+2a - 2b+c and then proceed to (u-a+b, v+a-b) . 
The number of such paths is, by the Reflection Principle, 



,u-a+b+v+a-b _ ,u+v x 
v+a-b- (2a-2b+c) v-a+b-o 



This technique may be extended to multiple reflections of any number 
of similar boundaries. 

2. The Number of Paths Touching a Boundary d Times 

Consider now a technique for counting the number of paths to 
(u,v) which touch boundary y = x + a, (a > 0) , d times, but do not 
cross it, and do not touch or cross boundary y = x + b, (b<0). This 
number is equal to the number of paths proceeding to (u,v) which touch 
y = x + a d times without crossing it less the number of paths touching 
y = x + a d times without crossing it that do touch or cross y = x + b. 
In [4] Milch and Waggoner developed a new technique called the Telescope 
Principle for counting the number of paths proceeding to (u,v) which 
touch y + x + a d times and do not cross it. This number is equal to 



•u+v-d, _ .u+v-d,. 
l u+a-r 1 u+a ; 



(5-4) 



To compute the number of paths to (u,v) which touch y = x + a 
d times, do not cross it, and touch or cross y = x + b, it is necessary 
to examine several disjoint subsets of that set of paths. 

If the number of paths in the set is N then 
N = + N 2 + N 3 where 

= the number of paths among these paths that touch 

or cross y = x + b before touching y = x + a 

= the number of paths among them that touch y = x + a 

d times without crossing it before touching or crossing 
y = x + b 
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( A 'n) 




> <0 
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Figure 



d-1 

N-. = Z (the number of paths among them that touch 
3 i=l 

y = x + a exactly i times and do not cross it before 
touching or crossing y = x + b) 

To find consider a path to (u,v) that crosses y = x + b 
first, then touches y = x + a d times without crossing it. Such 
a path is shown in Figure 5. If the path beyond its first point of 
contact with y = x + b is reflected about y = x + b, the reflected 
path (dotted line) will proceed to (v-b, u+b) and touch the reflection 
of y = x + a, i.e. , the line y = x + 2b-a, d times without ever 
crossing it. The Telescope Principle is then applied to the reflected 
path. Therefore the number of such paths is equal to 

[number of paths to (v-b, u+b-d) hitting neither y = x+a 
nor y = x+2b-a-d] 

- [number of paths to (v-b, u+b-d) hitting neither y = x+a 
nor y = x+2b - a-d+1] 

This number is obtained from (5-2) yielding the formula for N^: 

N 1 “ 1 u+a-b-1 ' u+a-b 1 " 1 ( u + 2a-b-l> 



,u+v-d \ -i , r /U+v-d v _ ,u+v-d ^ •. 
l u+2a-b 1 u-2a+3b l u-2a+3b+l J J 



(5-5) 



Next the number of paths, , proceeding to (u,v) which touch 
y = x + a d times but do not cross it then touch or cross y = x + b 
is counted. Such a path is shewn in Figure 6. If the Telescope 
Principle is applied to this path about boundary y = x + a, its 
equivalent path (dashed line) can be seen to meet the following 
criteria: 
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(A'n) (A'p-n) 
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Figure 



(i) it must touch but not cross line 

y = x + a + d-1 first and then must 
touch or cross line y = x + b + d. 
It may touch y = x + a + d-1 again 
after touching y = x + b. 

(ii) it must go to (u-d,v) . 

The number of paths so constrained is 



where 



n 2 = a- b- c + d 



r ,u+v-d « _ ,u+v-dx , 

1 v+a-b-1^ v+a-b^ J 

r ,u+v-d x _ ,u+v-d x , 
1 u+2a-b-l J u+2a-b 



(5-6) 



A = the number of paths to (u-d,v) which touch or 

cross y = x + a + d-1 first, then touch or cross 

y = x + b + d. From (5-1) , A = ( u ^" v ^ ,) . 

J v+a-b-1 

B = the number of paths to (u-d,v) which touch or 

cross y = x + a + d first, then touch y = x + b + d. 
Frcm (5-1) ,B = (^^) . 

C = the number of paths to (u-d,v) which touch or 
cross y = x + a + d-1 first, touch or cross 
y = x + b + d second, then touch or cross 
y = x + a + d. From (5-3) , C = £^_ d+1 ) • 

D = the number of paths to (u-d,v) which touch or 
cross y = x + a + d first, touch or cross 
y = x + b+d second, then touch or cross 
y = x + a+d again. D must be added since 
it has been subtracted in both B and C. 

Kan (5-3), D- C^b-d’- 
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Those paths which proceed to (u,v) and touch y = x + a (d-1) 
times before touching y = x + b, touch y = x + a a total of d times, 
and do not touch y = x + a + 1 are counted next. Figure 7 shews a 
path of this type. Applying the Telescope Principle to this path 
about boundary y = x + a its equivalent path (dashed line) can be 
seen to be constrained as follows: 

(i) it must touch boundary y = x + a + d-2 first, 
touch y = x + b + d-1 second, and then touch 



y = x + a + d-1. 

(ii) it must go to (u-d,v) . 

(iii) it may not touch y = x + a + d-1 until it has 

touched y = x + b + d-1. 

(iv) it may not touch boundary y = x + a + d. 

The number of such paths is equal to 

A' - B 1 - C' + D 1 = (“tT 3 . ,) - ( u ^- d . .) 

u+2a-b-2 u+2a-b-l 

, ,u+v-d > ,u+v-d v . 
u+2a-b-l' " u+2a-b J 



(5-7) 



where 



A' = the number of paths to (u-d,v) which must touch 

y = x+a+d-2 first, touch y = x+b+d-1 second, then 

touch y = x+a+d-1 . Fran (5-3) , A' = C-2a+b-d+2^ 

B' = the number of paths to (u-d,v) which must touch 

y = x+a+d-1 first, touch y = x+b+d-1 second, then 

touch y = x+a+d-1. From (5-3) , B' = ( U+ T , ) . 

v-2a+b-d+l 

C' = the number of paths to (u-d,v) which must touch 



y = x+a+d-2 first, touch y = x+b+d-1 second, then 
touch y = x+a+d. Fran (5-3), C' = ^-2a+b-d+l^ ’ 
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(A'n) (A'p-n) 
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Figure 7 



D' = the number of paths to (u-d,v) which must touch 
y=x+a+d-l first, touch y=x+b+d-l second, then 
touch y=x+a+d. D' must be added since it has 

been subtracted in both B' and C ' . From (5r3) , 

, ,u+v-d } 

V-2a+b-d ; * 

When the numbers of those paths proceeding to (u,v) which 
touch y=x+a 1,2, . . . > (d-2) times, respectively, before touching 
y=x+b, touch y=x+a a total of d times, and do not touch y=x+a+l 
are computed similarly, each is found to be equal to (5-7) . 
Therefore the total number of paths in is 

N 3 " " S - 1) u+2a-b-2 - 2 (d - 1) Wl 1 + (d - 1) <u + 2a-b> 



This may be written as 



N 3 " ^ >C^b-2» 



,u+v-d > , 
u+2a-b-l J 



, A ,x f ,u+v-d . ,u+v-d v, 
- W -11 u+2a-b-l - u+2a-b 



(5-8) 



Now N = + N^ + N^, and from (5-5) , (5-6) , and (5-8) 



m = f u+v-< 3 \ _ ,u+v-d> ,u+v-d v _ .u+v-d v 

u+a-b-1 u+a-b l u-2a+3b ; l u-2a+3b+l ; 



,u+v-d . ,u+v-d> M .,v r , u+v-d > ,u+v-d 

+ W-b-l* ' v+a-b ’ (d+1) u+2a-b-l ‘ UJ 1 



(5-9) 



/J T \ r ,u+v-d > ,u+v-d > , 
+ W’ 11 u+2a-b-2 - ut-2a-b-l 



Combining the results of (5-4) and (5-9) , the number of paths 
to (u,v) which touch boundary y=x+a d times, do not cross it, and do 
not touch boundary y=x+b is 
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r ,u+v-cL _ ,u+v-d, 1 r ,u+v-d 

u -‘- ^ ' u+a n l W-b-l 



v u+a-l 

,u+v-d 



) _ ( u+v-d }] 
‘ u+a-b' J 



,u+v-d 



,u+v-d 



-rr rv_ ^ ) - ( UTV ^ )] - [( UTV-U ) _ ( u+v-d ) i 
1 u-2a+3b' 'u-2a+3b+l ; J 1 'v+a-b-l ; 'v+a-b ; J 



, r ,u+v-d > ,u+v-d n M i > r /U+v-d . ,u+v-d , , 
+d u+2a-b-l " u+2a-b “ (d_1) [ ( u+2a-b-2 ) " ( u+2a-b-l )] 



,u+v-d 



u+v-d 



,u+v-d 



(5-10) 



C. THREE DIMENSIONAL COUNTING 

The number of paths from the origin to sane point (u,v,w) of the 
three supply ship model which do not touch any of the boundaries of 
Chapter III will be counted next. This number is equal to the total 
number of paths fran the origin to (u,v,w) , irrespective of boundaries, 
less the number of paths which touch at least one of the boundaries. 

The total number of paths fran the origin to (u,v,w) irrespective of 
any boundaries is the number of possible combinations of u steps in 
the x direction, v steps in the y direction and w steps in the z 
direction. This number is the multinomial coefficient, ( u+v+w ) = 

U /V 

( U+v+w ) _ l/u!viw! 

U/V/W 

In order to count the number of paths from (0,0,0) to (u,v,w) which 
contact boundary planes, it will be necessary to extend the Reflection 
Principle to three dimensions. By the Reflection Principle ( U "J‘ V+W ) / 

Vn~3i ,V ,U-a 

(U+v+w J and ( u+v+w ) are the number of paths that have at least 
v-a,u+a,w u,w-a,v+a ^ 

one common point with the planes x=z+a, y=x+a, and z=y+a, respectively, 

among the ( u+v+w ) paths connecting the origin and the point (u,v,w) . 

U /V ,w 

If N' is the number of paths from the origin to (u,v,w) which touch 
at least one boundary plane then 

N’ = E' + F’ + G’ - H' - I' - J' + K' 



where 
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E' = the number of paths touching x = z + M 

F' = the number of paths touching y = x + N 

G' = the number of paths touching z = y + P 

H' = the number of paths touching x = z + M and y = x + N 

I' = the number of paths touching x = z + M and z = y + P 

J' = the number of paths touching y = x + N and z = y + P 

K' = the nurber of paths touching x = z+M, y = x + N, and z 

By direct application of the Reflection Principle 

, = (U+V+W ) 

Wm,v,u-M 



= y + P. 



F' 

G' 



.U+V+W v 

v-N,u+N,w 



,u+v+w , 
u,w-P,v+P^ 



H ' is the sum of two numbers; the number of paths which strike x = z + M 

first and the number striking y = x + N first. I' and J' are the sum of 

similar terms. The formula for J' is explained below and H' and I' are 

listed in tabular form. The reflection of the plane y = x + N about the 

plane z = y + Pisx = z- (N + P). Any path proceeding to the point 

(u,v,w) , initially touching the plane z = y + P (one or more times) , and 

then proceeding to the plane y = x + N, is reflected beyond its first 

point of contact with z = y + P about that plane. The reflected path 

touches or crosses the reflection of y = x + N, i.e., the plane 

x = z - (N + P) , and proceeds to the point [ (u-N) , (v-P) , (w + N + P) ] . 

The number of such paths is ( u+ ^ +w ) . 

^ u-N, v-P 

The reflection of plane z = y + P about y = x + Nisx=z - (N + P). 

Any path proceeding to the point (u,v,w) , initially touching the plane 
y = x + N (one or more times) , and then proceeding to the plane z = y + P 
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is reflected beyond its first point of contact with y = x + N about 
that plane. The reflected path touches or crosses the reflection of 
z = y + P, i.e. , the plane x = z - (N + P) , and proceeds to the point 
Hu - (N + P» , (v + N) , <w + P)] . Therefore J' = <£££*> + C^. p > • 
H' and I' are similarly computed and are tabulated below. 



Initial Plane 
Second Plane 
Reflected Plane 
Reflected Point 
Number of Paths 



Table for H' 
y=x+N 
x+z+M 
z=y- (M+N) 

(u-N , v+MfN ,w-M] 

.u+v+w . 
u-N,w-M 



X=z+M 
y=xtN 
z=y— (M+N) 

[u+M,V+N,w- (M+N] 

,u+v+w . 
U+M,v+N 



, .u+v+w . ,u+v+w . 
V u-N,w-M ; V u+M,v+N ; 



Initial Plane 
Second Plane 
Reflected Plane 
Reflected Point 
Number of Paths 



Table for I' 

z=y+P 
x=z+M 
y=x- (M+-P) 

[u+M,v- (MtP) ,v+P] 

.u+v+w , 
u+M,w+P 



.u+v+w . .u+v+w 

u+M,w+P v-P ,w-M' 



x=z+M 
z=y+P 
y=x- (M+P) 

[u+M+P , v-P ,w-M] 

.u+v+w . 
v-P ,w-N 



K may be ccmputed in similar fashion by permiting all three boundary 
planes and applying the Reflection Principle twice after each permutation. 
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